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We discuss whether an appropriately denned dimensionless scalar function might be an acceptable 
candidate for the gravitational entropy, by explicitly considering Szekeres and Bianchi type Vlh 
models that admit an isotropic singularity. We also briefly discuss other possible gravitational 
entropy functions, including an appropriate measure of the velocity dependent Bel-Robinson tensor. 

I. INTRODUCTION 

Penrose [1] has argued that the initial cosmological singularity must be one of low entropy in order to explain 
the high isotropy of the observed universe and to be consistent with the second law of thermodynamics. Since the 
matter was presumably in thermal equilibrium at the initial singularity, this implies low entropy in the gravitational 
field. Penrose has also conjectured that such a gravitational entropy should be related to a suitable measure of the 
Weyl curvature. The search for a suitable candidate for the "gravitational entropy" is therefore of current interest, 
particularly in the approach to the initial cosmological singularity [2,3]. In the quiescent cosmological paradigm [4] 
the universe began in a highly regular state and subsequently evolved towards irregularity. 

Penrose [1] originally proposed that the the Weyl tensor is zero at the big bang singularity, implying the subsequent 
evolution is close to a Friedman-Robertson- Walker (FRW) model. However, this requirement is too strong. For 
example, in perfect fluid spacetimes Anguige and Tod [5] have proven uniqueness results that show that if the Weyl 
tensor is zero at the big bang, then the spacetime geometry must be exactly FRW in a neighbourhood of the big- 
bang. This has motivated the idea that some appropriate dimensionless scalar is asymptotically zero. The search for a 
gravitational entropy then reduces to a search for this scalar function. Quiescent cosmology and the ideas of Penrose 
provide the motivation of the definition of an isotropic singularity [6]; essentially a spacetime admits an isotropic 
singularity if the 'physical' spacetime is conformally related to an 'unphysical' spacetime such that there exists a 
time function T with the property that at T = the conformal factor vanishes (corresponding to the cosmological 
singularity) but that the conformally related metric is regular. It was proven in [6] that in the class of models with 
an isotropic singularity P — > as T — > 0, where P is the ratio of the Weyl curvature squared to the Ricci curvature 
squared (see Eqn. (3) below). 

Recently, Lake and Pelavas (LP) [7] considered a class of "gravitational epoch" functions that are a dimensionless 
scalar field constructed from the Riemann tensor and its covariant derivatives only (such as, for example, P). They 
discussed whether such functions can act as a "gravitational entropy" by determining whether it is monotone along a 
suitable set of (smooth) timelike trajectories. In particular, LP considered the set of homothetic trajectories of a self- 
similar spacetime (since such spacetimes are believed to play an important role in describing the asymptotic properties 
of more general models [8] ) . They showed that the Lie derivative of any " dimensionless" scalar along a homothetic 
vector field (HVF) is zero, and concluded that such functions are not acceptable candidates for the gravitational 
entropy. They suggested considering other options for a " gravitational epoch" function. Other dimensionless scalars 
constructed from the Riemann tensor and its covariant derivatives have been considered (see, for example, [3]). Other 
alternatives, such as those involving the Bel- Robinson tensor, were suggested in [9]. 

In this paper, by explicitly considering classes of Szekeres and Bianchi VI^ models that admit an isotropic singularity, 
we revisit the conclusion of LP that P (for example) is not an acceptable candidate for the gravitational entropy. 
First, we take the view, unlike that taken in LP, that a purely gravitational entropy selects spacetimes as being of 
cosmological interest according to a thermodynamic principle. For example, in General relativity (GR) there exist 
solutions for which the physical energy density is negative; however, we do not disregard GR or the notion of energy 
density within GR - rather, we use the criterion of negative energy density to characterize those solutions which 
are not physical. Second, homothetically self-similar spacetimes represent asymptotic equilibrium states (since they 
describe the asymptotic properties of more general models [8]), and the LP result is perhaps consistent with this 
interpretation since the entropy does not change in these equilibrium models, and perhaps consequently supports 
the idea that P (for example) represents a "gravitational entropy". Therefore, we will investigate the behavior of 
P asymptotically as the self-similar solution is approached in cosmological models. We also consider various other 
options for a purely "gravitational entropy" , including an appropriate measure of the Bel- Robinson tensor. 



1 



II. SZEKERES MODEL 



We consider the class II Szekeres solutions, which are spatially inhomogcncous models with irrotational dust as 
a source. A comprehensive analysis of the Szekeres models can be found in [10]. In the notation of Goode and 
Wainwright we set k — = /3_ , the line-element in comoving coordinates has the form 

ds 2 = T 4 [-dT 2 + dx 2 + dy 2 + (A - p + T 2 fdz 2 ] (1) 

where 

A = a{z) + b(z)x + c(z)y - 5(3+(z)(x 2 + y 2 ). (2) 

The fluid 4-velocity is u = T~ 2 -j^ and the energy density is [i = 12T~ 6 [1 — (f3 + / A)T 2 ]~ X . These cosmological models 
admit an isotropic singularity at T = 0; moreover, it has been shown [11] that the general Szekeres class with /3_ = 
also admits an isotropic singularity. If f3 + = in (1), then we obtain the associated flat FRW dust solution. 
The standard gravitational epoch function, P, for (1) is 

C abcd C abcd _ AT^X 
R ab R ab 3 A 2 ' ( ' 

This has been shown [12] to behave appropriately in these models, i.e. P is monotonically increasing away from 
the isotropic singularity. Perhaps an alternative choice is to use the Bel-Robinson tensor [9] to construct a velocity 
dependent gravitational epoch function. Using the fluid 4-velocity, we construct the positive scalar 

24/3 2 

W = T abcd u*uW = T&{p+T2 + _ A)2 (4) 

which has the same units as C abc dC abcd , and hence to obtain a dimensionless scalar we normalize by the square of 
^ = T a bu a u b , i.e. 



W _ T 4 l3l 
/i 2 ~ 6A 2 



P=^ = ^T- (5) 



Noting that in this case we have P — P/8 then P also behaves appropriately as the isotropic singularity is 
approached. This relationship is a consequence of the following two facts. First, the magnetic part of Weyl, H ab = 0, 
therefore C abc( iC abcd is equivalent to W modulo a positive constant. Second, the Ricci invariant R ab R ab = R 2 
and for dust R — \i. In these models, the choices are limited for constructing dimensionless ratios of zeroth order 
invariants since all Carminati-McLenaghan (CM) invariants 1 can be expressed in terms of the Ricci scalar R and 
Re(wi) - C abcd C abcd , 

R = — = 2AP + ffil 

T 6 (A-/3 + T 2 )' Wl T*(A-f3+T 2 ) 2 [ ' 

(7) 

with syzygies for the Ricci invariants, 

n = ±R 2 , r 2 = l A R\ r 3 = T i I i 1 ' 4 . (8) 
Since these spacetimes are Petrov type D, the Weyl invariants satisfy 6w 2 = ivf . The mixed invariants give syzygies 
mi = 0, m 2 = jgWiR 2 = m 3 , m 4 = — j^wiR 3 , m 5 ~j^w 2 R 2 - (9) 

To study the Lie derivative of P (or P) along a HVF we first recall a well known result [8] . Any FRW model with an 
equation of state p = — admits a timelike HVF. In addition, only the flat FRW models with an equation of state 
p = a[i and power law dependence on the scale function admit a HVF. If /3+ is set to zero in (1), then the resulting 



1 We note that in this case any complex CM invariants have vanishing imaginary part. 
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£ • £ = ( f) n-T 2 + ? + [!- £ /3+(i? 2 + T 2 )] 2 z 2 } (12) 



flat FRW model is in non-standard coordinates; although there exists a HVF, finding it in this coordinate system is 
difficult due to the presence of the functions a, b and c. Determining a coordinate transformation into more familiar 
coordinates, where HVF's can easily be found, is also difficult. We use the following simplifying assumptions in (1), 
a = 1, b = 0, c = 0, and consider deviations from flat FRW by redefining /3 + — > e(3 + for small e. It can be shown that 

, <t> ( m d d d d \ 

Z=-3{ T df +X d- X +y d-y +Z d- Z ) (10) 

is a HVF in the associated fiat FRW model and deviates 2 from homotheticity in (1) to first order in e. Since 
limT^o(I j {(7o6 — 20<7 a b) = we find that (10) becomes a HVF as the isotropic singularity is approached. In this limit 
then L^P = 0, as expected [7], and so P approaches a constant. This type of behavior is desirable if we would like 
to interpret P, in some sense, as a gravitational entropy possessing a critical point at early times in the evolution. 
Moreover, requiring that P be monotonically increasing at early times along a timelike £ (I^P > 0) places restrictions 
on j3+. Since 

and 

Zj 5 

where R 2 = 5(x 2 + y 2 ), then as T — > 0, £ remains timelike if R = and z — > subject to the requirement that 
lim^o zj3 + be bounded. Assuming z(3 + is analytic near z = gives the form 

/?+ = - + 6i + 6 2 z + • • • (13) 

z 

To leading order in e we factor (3 + from (11) and use (13) to obtain lim z ^o(2 + z/3 + /[3 + ) = 1 thus L^P > along 
timclike £. 

III. BIANCHI VI H MODEL 

It has been shown [13] for the Bianchi VI^ class that a choice of parameters can result in the quasi-isotropic stage 
beginning at the initial singularity giving rise to an isotropic singularity for these spacetimes. In the notation of [13], 
we set a s — and a m = 1, so that the line-element in conformal time coordinates is 

ds 2 = r 4 /( 3 ^ 2 > (-A 2 ^dT 2 + A 2 ^dx 2 + A 2 ^e 2r ^+^- 2 ^dy 2 + A^e^-^-^dz 2 ) (14) 

where 

^ - 1 + « - 1 - ^ «. - ^ «• - <* + ™ - * - 4(2 p ;^v (») 

These spacetimes have a perfect fluid source with equation of state p = (7 — l)/i, 1 < 7 < 2. The fluid 4-velocity is 
u = A^t- 2 '^- 2 ^ and the energy density is \i = l2A-~t(3j-2)- 2 T- 6 ~t/<- 3 ~>- 2 \ Since a s has been set to zero, then 
the isotropic singularity occurs at t = 0. The parameter a c determines the curvature of the spacelike hypersurfaces 
orthogonal to u, if a c = we obtain the flat FRW solution. We shall consider deviations about this flat FRW model 
by assuming a c is small. 

To leading order in a c , the gravitational epoch function for (14) is 



As a special case, if (3+ = C/z 2 then (10) is also a HVF of (1); therefore, as a consequence of [7], L^P = 0. 
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which is positive and P — > as r — ► + . Using the Bel- Robinson tensor and the energy density we find 

~ 7 2 (37- 2) 2 t 4 2 3 ^ M7 x 

P= 6(7-2)2 Q c + Q(«c)- (17) 

Consequently to leading order in a c we have that P = [3(7— l) 2 + l]P/8, and again P and P are directly proportional. 
Unlike the Szekeres class above, the magnetic part of the Weyl tensor with respect to u does not vanish here unless 3 
7 = 4/3. This relationship between P and P becomes evident if we consider the expansions for the relevant invariants 
of the electric and magnetic parts of Weyl 

F F ab_ 24V „? + n(„*\ H H« b - 2 7 2 (3 7 + 2)(3 7 - 4) 2 3 4 

EabE ~ (7 _ 2) 2 ( 3 7 _ 2) 2 T 8/(3 7 -2) a C + ( a c), ^a 6 P - (JJ _ 7 )3(3 7 _ 2 )2 T 6(2- 7 )/(3 7 -2) ^ + °^c)- ( 18 ) 

For a c small, H ab H ab ~ and C afccd C af>c<i - E ab E ab ~ W, additionally Einstein's equations give R ab R ab = [3(7 - 
l) 2 + l]/x 2 , the relationship now follows. 
As a c — > or if r — > + , the vector 

, (3 7 - 2) / d d d d \ 

gives (L^gab — 2<p9ab) — * separately in both limits; therefore, ^ is a HVF in the flat FRW limit or as the isotropic 
singularity is approached. Setting i = i/ = z = 0we obtain for the magnitude of £ 



£ . f = _ ^ 2 fc_2l! T 67/(37-2) A 2( 7 -l) . 



(20) 



thus by continuity ^ will be timelike in some neighborhood of x = y = z = arbitrarily close to the isotropic 
singularity. To leading order in a c , the behavior of P along £ is given by 

which is positive for 1 < 7 < 2; hence P will be monotonically increasing at early times along the timelike £, and as 
t — > + then P will tend to a constant (which in these models is zero). 

We now show that the Weyl curvature hypothesis does not necessarily put restrictions on the Petrov type. In this 
class of spacetimes the Weyl invariants of CM have vanishing imaginary parts, and in general do not always satisfy 
the syzygy wf — 6w| = 0; therefore, the metric(14) is almost always Petrov type I. However for 4/3 < 7 < 2 there 
always exists a time r* where the Petrov type specializes to either II or D; this is given by 

r. = I . 2 ~ 7 , (22) 
2^(37-4)' 

otherwise the Petrov type is I. When 7 = 4/3 the syzygy is satisfied and by choosing the aligned Newman-Penrose 
tetrad 



we find that this is in fact Petrov type D for all r > 0. It would appear that an intermediate algebraic specialization 
of the Weyl tensor during the evolution does not affect the increasing anisotropy that is indicated by the gravitational 
epoch function (16); indeed, these spacetimes begin with small anisotropy close to the isotropic singularity and 
approach the anisotropic vacuum plane wave metrics at late times. 

As is shown in [7], any dimensionless ratio of invariants is constant along a HVF. Therefore, depending on its 
monotonicity, it may also serve as a gravitational epoch function when the cosmological model admits an asymptotic 



3 H a b will also vanish if a c = 0. 
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timclike HVF. Here we illustrate this point by considering a dimensionless ratio of differential invariants; to second 
order in a c we have 

D _ VaC bcde \7*C bcde _ 40 (3 7 - 2) 2 r 4 2 3 

n_ V a P 6c V«P> " 9 ( 7 -2)2[9( 7 -l)2 + 5 ] ac + U(acj (24j 

and 

Clearly as r — > + , Pi — » and from (25) it is also monotonically increasing for 1 < 7 < 2. Nevertheless, the invariants 
of (24) diverge in a similar manner to the invariants of (16), i.e. as r — > + , V a C{, c( zeV a C bc<ie and V a i?b c V a i? bc — > —00. 



IV. DISCUSSION 



We have considered a class of spatially inhomogeneous Szekeres solutions with the line-element (1). We show that 
there exists a HVF (10) as the isotropic singularity is approached, and show that L^P = in this limit, so that P 
approaches a constant as expected [7]. Moreover, to leading order we show that L^P > along an approach to the 
singularity in which the HVF £ remains timelike; i.e., P is monotonically increasing at early times along £. We then 
considered a class of Bianchi VI/j models with an isotropic singularity with line-element (14), parameterized by a c 
(which measures deviations about the flat FRW model). Assuming a c is small, to leading order we display a vector 
£ which is a HVF in the flat FRW limit, and show that the gravitational epoch function P — > as the isotropic 
singularity is approached and that P is monotonically increasing at early times along timelike £. 

Therefore, in the isotropic singularity cosmological models we have studied we have found that P — > asymptotically 
as the self-similar cosmological model is approached, in support of the idea that these homothetically self-similar 
spacetimes represent asymptotic equilibrium states. Moreover, we have provided evidence that P is monotonically 
increasing as the models evolve away from these equilibrium states, which perhaps lends support to the idea that P 
represents a " gravitational entropy" . 

We also found that for both the Szekeres models and the Bianchi VI h models (to leading order in a c ), the standard 
gravitational epoch function P, and the normalized Bel-Robinson epoch function P, are proportional (and hence 
equivalent as gravitational epoch functions). The question remains as to whether this will be true for all models 
with an isotropic singularity. In general, for a perfect fluid source we have that P - (P 2 + H 2 )/^ 2 and P ~ 
(E 2 - H 2 ) / (fi 2 + 3p 2 ) . Assuming an equation of state of the form p = ay gives P ~ (P 2 - Pf 2 )/[(1 + 3a 2 )/z 2 ]. Clearly 
if H 2 is negligible with respect to P 2 , then P and P will be effectively proportional. It may also be of interest to 
consider cosmological models where P and P differ. In particular, whenever the Petrov type is III, N or O then all 
zeroth order Wcyl invariants vanish, and hence P vanishes but P does not necessarily vanish. An example of such 
cosmological models are the Oleson [14] solutions, which are Petrov type N with a perfect fluid source. In these 
models P vanishes but P does not; it is of interest to determine if these models can admit an isotropic singularity. 

A classic problem in cosmology is finding a way to explain the very high degree of isotropy observed in the cosmic 
microwave background. In GR cosmological models admitting an isotropic singularity are of zero measure, so that 
isotropy is a special rather than generic feature of cosmological models. Hence, a dynamical mechanism which is able 
to produce isotropy, such as inflation, is needed. However, inflation requires sufficiently homogeneous initial data in 
order to begin [15]; hence the isotropy problem remains open to debate in standard cosmology. Recently, it has been 
argued that an isotropic singularity is typical in brane world cosmological models [16]. Hence brane cosmology would 
have the very attractive feature that it provides for the necessary sufficiently smooth initial conditions which might, 
in turn, be consistent with entropy arguments and the second law of thermodynamics. 



V. ACKNOWLEDGEMENTS 



We would like to thank K. Lake for helpful discussions. AAC is supported by NSERC of Canada. 



5 



[1] R. Penrose, in General Relativity, an Einstein centenary survey, eds. S.W. Hawking and W. Israel, Cambridge Univ. Press 
(1979). 

[2] P. Tod, Class. Quant. Grav. 7, L13 (1990). 

[3] S. W. Goode, A. A. Coley and J. Wainwright, Class. Quant. Grav. 9, 445 (1992). 

[4] J. D. Barrow, Nature 272, 211 (1978). 

[5] K. Anguige and K. P. Tod, Annals Phys. 276, 257 (1999). 

[6] S. W. Goode and J. Wainwright, Class. Quantum Grav. 2, 99 (1985). 

[7] N. Pelavas and K. Lake, Phys. Rev. D. 62, 044009 (2000). 

[8] B. J. Carr and A. A. Coley, Class. Quant. Grav. 16, R13 (1999) and references therein. 

[9] N. Pelavas and K. Lake, gr-qc/9811085 version 1. 
[10] S.W. Goode and J. Wainwright, Phys. Rev. D. 26, 3315 (1982). 
[11] S.W. Goode and J. Wainwright, Mon. Not. R. Astron. Soc. 198, 83 (1982). 
[12] S.W. Goode, A. A. Coley and J. Waingwright, Class. Quant. Grav. 9, 445 (1992). 
[13] J. Wainwright and P.J. Anderson, Gen. Rel. Grav. 16, No. 7, 609 (1984). 

[14] H. Stephani,D. Kramer, M. MacCallum and E. Herlt Exact Solutions to Einstein's Field Equations, Second Edition, Cam- 
bridge Univ. Press (2003). 
[15] E. W. Kolb and M. S. Turner, The Early Universe, (Addison - Wesley, 1990). 

[16] A. A. Coley, Phys. Rev. D 66, 023512 (2002); A. A. Coley, Class. Quant. Grav. 19, L45 (2002); A. A. Coley, Y. He and 
W. C. Lim, Class. Quant. Grav. 21, 1311 (2004). 



G 



